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ABSTRACT 

We measure the radial profile of the CO (1-0) to H2 conversion factor (Xco) in NGC 628. The Ha 
emission from the VENGA integral field spectroscopy is used to map the star formation rate surface 
density (Y]sfr)- We estimate the molecular gas surface density (S//2) from ^sfr by inverting the 
molecular star formation law (SFL), and compare it to the CO intensity to measure Xco- We study 
the impact of systematic uncertainties by changing the slope of the SFL, using different SFR tracers 
(Ha vs. far-UV plus 24^m), and CO maps from different telescopes (single-dish and interferometers). 
The observed Xco profile is robust against these systematics, drops by a factor of 2 from R ^ 7 kpc to 
the center of the galaxy, and is well fit by a gradient Alog(Xco) = 0.06 ± 0.02 dex kpc~^. We study 
how changes in Xco follow changes in metallicity, gas density, and ionization parameter. Theoretical 
models show that the gradient in Xco can be explained by a combination of decreasing metallicity, 
and decreasing with radius. Photoelectric heating from the local UV radiation field appears to 
contribute to the decrease of Xco in higher density regions. Our results show that galactic environ- 
ment plays an important role at setting the physical conditions in star forming regions, in particular 
the chemistry of carbon in molecular complexes, and the radiative transfer of CO emission. We cau- 
tion against adopting a single Xco value when large changes in gas surface density or metallicity are 
present. 

Subject headings: galaxies: ISM 



1. INTRODUCTION 

Measuring and studying molecular gas in galaxies is 
fundamental to understand star formation, and the phys- 
ical processes setting the balance between the differ- 
ent phases of the interstellar medium (ISM). Molecu- 
lar hydrogen (H2) amounts for the bulk of the mass 
in molecules in the universe, but its observable transi- 
tions are rarely excited at the typically cold tempera- 
tures ('--'lO K) of the gas inside giant molecular clouds 
(GMCs). To overcome this observational difficulty, the 
second most abundant molecule in GMCs, the carbon 
monoxide molecule ^^C^^O (hereafter CO), is typically 
used as a proxy for estimating the total mass in H2. 
Rotational CO transitions, observed at millimeter (mm) 
wavelengths, can be easily excited under the typical den- 
sity and temperature conditions in GMCs, and are there- 
fore bright enough to be detectable in single molecular 
clouds in the Milky Way (MW), galaxies within the Lo- 
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cal Group, across the disks of nearby galaxies out to dis- 
tances of ~ 10 Mpc, and even out to high redshifts by 
integrating the emission over whole galaxies. 

Using CO emission to estimate the H2 mass requires 
knowledge of the CO intensity to H2 column density con- 
version factor: 



X: 
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I{CO) ' 



(1) 



which for Galactic molecular clouds in the vicinity of the 
Sun, has typical values of 2 — 4 x 10^° cm~^(K km s~^)~^ 
(|Kennicutt fc Evansl 120121 and references therein). Or 
alternatively the CO luminosity to total mass conversion 
factor aco = Mh2/ L{CO) which unlike Xco includes a 
factor of 1.36 for the contribution of Helium to the total 
masjn. 

Constraining the appropriate value of this conversion 
factor and establishing how it changes under the different 
physical conditions present across different environments 
inside and across galaxies is of paramount importance 
if astrophysical interpretations regarding the molecular 
ISM are to be drawn from CO data. Of particular inter- 
est is the metallicity dependance of Xco, as deep obser- 
vations with current facilities, and the advent of new and 
more powerful telescopes like ALMA and CCAT, allow 
the detection of CO in dwarf galaxies, the outskirts of 
the disks of massive spirals, and high redshift systems, 
where the heavy element abundance is expected to be 
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low. 

A series of studies have used different techniques to 
explore the possibility of a changing Xqo across dif- 
ferent types of galaxies in the nearby universe. These 
include virial mass measurements of individual GMCs 
in t he Milky Way , the local group , and nearby spirals 
(e.g. I Wilson! IT995I: IBlitz et al.l ]2007l: iBolatto et al.l l2008l: 
iFukui fc Kawamurall2010l and references therein), esti- 
mating the molecular gas mass from dust far-IR emis- 
sion modeling while constraining the dust-to-gas ratio 
and the contribut i on fro m atomic hydrogen (Israel 1997; 
ILerovet all IMIIMI . and using the star formation 
rate (SFR) under the assumption of a known molecu- 
lar gas depletion timescale to estimate the a mount of H2 
(jSchruba et al.l [20I2I: iMcQuinn et al.l IMl) . A consis- 
tent picture seems to arise from these studies, in which 
Xco shows higher values for lower metallicity systems. 
The difference can be dramatic for the lowest metallic- 
ity dwarfs in the local universe, where the conversion 
factor can be 10 to 100 times higher than in the Milky 
Way. This increase is most likely driven not only by 
a decrease in the carbon and oxygen abundances, but 
mainly by a drop in the optical depth within GMCs 
due to a lower abundance of dust. The latter trans- 
lates in the CO/C^ dissociation boundary moving in- 
wards within these clouds, leaving behind large envelo pes 
of "CO dark" molecular gas fe.g. IBolatto et al1ll999[ ). 

On the other hand, studies of molecular gas in merg- 
ing and starburst galaxies (typically ultra-luminous in- 
frared galaxies, ULIRGS) based on virial mass mea- 
surements, dust emission modeling, and column den- 
sity estimation from optically thin transitions of CO 
isotopes, find Xco values which are factors of a few 
lower than the typical MW values in the solar v icinity 
JWild ct al. 19 92; Shi er ct al. 19 94: Mauersberger~eral] 
1996; Solomo n et all Il 997: D ownes fc SolomonI 119981 
Bryant fc Scovilld 11999c ,Meier et al.l I2010D . The 
same effect is observed for ULIRGS and sub- 
mn i galaxies (SMGs) at h igh r edshi ft (z > 1) 
bv iSolomon fc Vanden Bou^ (|2005D and iTacconi et al.l 
(poos). This effect is thought to be caused by the im- 
pact of higher gas temperatures and stronger turbulence 
on the brightness temperature of the CO line and the 
escape probability of CO (1-0) photons. Since the C0(1- 
0) transition is typically optically thick, the broadening 
of the line-width (Av) induced by higher levels of tur- 
bulence in these high densit y environments pro motes the 
escape of CO(l-O) photons (jShettv et al.ll2011l ). Interest- 
ingly, the Xco dependance with metallicity mentioned in 
the last paragraph has also been observed in a sample of 
more "norma l" star forming galaxies at high redshift by 
iGenzel et all (|2012) using a method similar to the one 
used in this work. 

These observational efforts to measure changes in 
Xco across different environments have been accom- 
panied in the last few years by detailed theoreti- 
cal modeling attempts to understand how CO ra- 
diative transfer depends on the physical conditions 
of the ISM. A series of studies using analytic mod- 
els, numerical simulations, and combinations of both, 
have examined the dependance of Xco with metal- 
licity, gas tempera ture, gas dynaruics, and the local 
radiation field (e.g. [ Kr umholz et al.' 2011; 'Shett vet al.l 
[2011;, Narayanan et al.F2012; .Feldniann ct al...2012al) . al- 



though see early work by IDickman et al.l ()1986l ) and 
IMalonev fc BlacH ([1988). 

Considering the fact that Xco changes from galaxy 
to galaxy depending on the average physical conditions 
of the ISM, it would not be surprising if it also changes 
within galaxies, depending on the local physical condi- 
tions present in different environments inside an individ- 
ual system. For years, this has been known to be the case 
in our own galaxy. The value of Xco has been shown to 
change as a function of galactocentric radius in the MW 
using a scries of diffe rent t echniques: dust emission mod- 
eling fSodroski et al.lll9"95l ). measurements of gamma-ray 
emiss i vity from cosmic-r a y gas interactions ()Digel et al.l 
[19961 IStrong et al.l [20M: lAbdoeFaD [Moh , and direct 
virial mass measurements of GMCs (Arim oto "eFalllTggel: 
lOka et all Il998). All these studies find a decrease in 
Xco from the typical values measured near the solar ra- 
dius (2 — 4 X 10^° cm~^(K km s~^)~^) towards smaller 
galactocentric radii, with the conversion factor reaching 
typical starburst /merger type values of 0.1 — 0.5 x 10^° 
cm~"'^(K km s^^)^^ in the Galactic Center. It is not 
clear from current studies if the Xco radial profile in the 
Milky Way follows a smooth gradient or if it is fairly con- 
stant across the Galactic disk and falls sharply at some 
intermediate radius. 

Measurements of the spatial distribution of Xco 
within spiral galaxies are scarce in the literature, and 
the subject remains highly unexplored from an obser- 
vational perspective. Observations of the optically thin 
CO isotopes ^■^CO and C^^O, and dust continuum emis- 
sion, have been used to show that Xco is a factor of 2-4 
lower than the typical solar vicinity MW values in the 
central region s of the nearby barred spirals NGC 694 6 
and Maffei 2 (|Meier fc Turned 120041: iMeier et al.l 12008V 
confirming the results from an earlier study of NG C 253, 
IC 342, Maffei 2, and NGC 6 9 46 by IWailetaLl {1993). 
Similar work by iVila-Vilarol (|2008| ) has shown a de- 
creased Xco factor in the central kpc of NGC 5194 
(a.k.a. M51a). Using an independent method based on 
estimating the dust mass distribution by conducting ra- 
diative transfer modeling of optical and near-IR images 
of galaxies, a nd assuming a constant dust-to-gas ratio, 
iRegaiJ (j2000l ) also finds low Xco values in the central re- 
gions (< 1.5 kpc) of NGC 1068, NGC 1530, NGC 2903, 
and NGC 6946. These results are consistent with what 
is observed in the MW. The radial profile of Xco out to 
large radii in spiral galaxies other than the MW remains 
largely unexplored. 

Only three published measurements of the Xco radial 
profile for spiral galaxies other than the Milky Way are 
known to the authors. First, the work of lArimoto et al.l 
(1996?) in NGC 5194 (M51a), who use CMC virial mass 
estimates to find that Xco follows a linear gradient inside 
one effective radius in t his galaxy. Their m easurements 
are based on the data of lAdler et al.l (jl992n which has a 
^^350 pc beam-size, and we do not consider virial masses 
measured on such physical scales to be reliable. Sec- 
ond, is the high resolu tion (beam-size ^ 20 pc) study of 
GMCs in M33 bv iRosolowskv et all ()2003D , who use virial 
masses derived from CO line-widths to find a flat Xco 
distribution across the galaxy, even in the presence of a 
0.8 dex change in metallicity across their sample. And, 
finally, the soon to be published dust modeling study 
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of a sample of nearby spirals by Karin Sandstrom (pri- 
vate communicatio n), which is discussed in §4 (see also 
lAniano et al.|[20T^ . The spatial distribution of the con- 
version factor in external galaxies is the subject of this 
paper, in which we study the radial profile of the Xqq 
across the disk of the nearby face-on Sc galaxy NGC 628. 

Our method to measure Xqo is based on the cor- 
relation observed over many orders of magnitude be- 
tween the SFR surface density and the surface den- 
sity of molecular gas in the ISM of star forming galax- 
ies. This relation, typically known as the "star for- 
mation l a w" (SFL) or the "Sch imdt-Kennicutt law" 
(jSchmidtl Il959t iKennicuttI I1998D . has been recently 
constrained in a spatially resol ved manner across the 
disks of nearby spiral galaxies (Kennic utt et al. 
lifiicl ct al. 2008; Blanc ct al. 2009; Vc rlev et al 



2007 



2010; 



2011 



Onodera ct al. 201Ck lSchruba et al.ll201lHLiu et al 
Rahman e t aL _2Q12(). Although debate persists regard 



ing the actual value of the SFL slope, which is subject 
to a series of systematic uncertainties related to back- 
ground subtractio n, cloud sanipling and fitting methods 
(see discussions in'Blanc et al.|[2009t IRahman et al.ll2012l ; 
ICalzctti ct al. 2012), its normalization is consistent with 
a depletion timescale for molecular gas of ~ 2 Gyr at the 
typical molecular gas surface densities observed across 
the disks of nearby spirals (Eg2 = 10 — 1 00 Mopc~^, 
e.g. iLerov et all 120081; IRahman et"aIll2Q12D . This rela- 
tion can be used to derive the molecular gas surface den- 
sity from the observed SFR surface density across the 
disk of a galaxy, therefore permitting the measurement 
of Xco by comparison to CO intensity maps. 

We describe the multi- wavelength datasets used in this 
work in §2, and our method for estimating Xco in §3. 
We test the robustness of our method to systematic un- 
certainties, by changing the assumed value for the slope 
of the molecular SFL, using different SFR indicators, and 
CO maps from different telescopes (both single-dish and 
interferometers). Our results are presented in §4, where 
we present the observed Xco radial profile of NGC 628 
and compare it to that derived from dust emission mod- 
eling using Spitzer+Herschel photometry (Karin Sand- 
strom private communication). In §5 and §6 we discuss 
the possible physical origins for the observed spatial dis- 
tribution in Xco J including the roles of metallicity, gas 
surface density, and the local UV radiation field. Finally 
we provide our conclusions in §7. For NGC 628, we as- 
sume an inclination oi i = 8.7° and a distance of 8.6 Mpc 
(jHerrmann et al.ll2008l ). 

2. DATA 

2.1. VENGA Integral Field Spectroscopy 

All the measurements of nebular emission lines used 
to estimate the metallicity, ionization parameter, dust 
extinction through the B aimer decrement, and the Ha 
SFR, are made on the Mitchell Spectrograph (formerly 
VIRUS-P) IFU data-cube of NGC 628 produced by the 
VENGA survey (jBlanc et al.ll2010( ). The data-cube sam- 
ples a rectangular area of 5.2' x 1.7' centered in the 
nucleus of the galaxy, and has a spatial resolution of 
5.6" full- width half-max (FWHM). The spectra covers 
the 3555A-6790A wavelength range with an instrumen- 
tal spectral resolution of ~ 5A FWHM. At the assumed 
distance of NGC 628 the VENGA spatial resolution cor- 



responds to ^235 pc, a few times larger tha n the typical 
sizes of large giant molecular clouds (< 60 pc lHever et al.l 
I2Q01D and indi vidual HII regions ion ized by single clus- 
ters (< 120 pc. iGutierrez et aDl2011[ ). The observations, 
data reduction, and calibrations will be described in an 
upcoming publication. 

2.2. GALEX far-UV and Spitzer MIPS 24 iini Data 

To ensure that our results are robust against the choice 
of the SFR tracer used, along with the VENGA Ha de- 
rived SFRs, we use a linear combination of far-ultraviolet 
(FUV) and mid-infrared (mid-IR) 24/xm emission. The 
former traces unobscured star formation, while the latter 
recovers dust obscured star formation by means of repro- 
cessed UV radiation reemitted as thermal IR emission 
from heated interstellar dust grains. A number of cali- 
brations have been proposed to use linear combinations 
of obscured and unob scured tracers to estimate the SFR 
(Calzctti ct al. 20 071 ILerov et al. 2008; Kennicut t et al 
2009t iHao et al.ll2011l: iMurphv et al.ll2011l; iLerov et al 

mm. 



We use the GALEX FUV ima ge of NGC 628 from the 
GALEX Nearby Galaxy Survev (|Gil de Paz et al.l[2?)0l 
which samples the 1350-1750 A wavelength range, has a 
point-spread-function (PSF) FWHM of 4.5", and is deep 
enough for us to measure the FUV flux at high signal-to- 
noise over the whole area of interest. To trace the IR dust 
emission we use the MIPS 24//m data of NGC 628, taken 
as part of the Spitzer Infrared Nearby Galaxy Survey 
(SINGS, iKennicutt et al.l I2003D and t he Local Volume 
Legacy survey (LVL. IDale erahl 120091 ). The map has a 
PSF FWHM of 6" and like the FUV map, it is deep 
enough for us to measure the 24/im flux over the whole 
area of interest in this paper. The data processing of 
both the FUV and 24/zm maps, including the suppression 
of backgro unds and th e mask ing of foreground stars is 
detailed in ILerov et al.l ()2012| ). 

2.3. HERAGLES, BIMA-SONG and GARMA CO Data 

We use three independent datasets to measure I{CO) 
across the disk of NGC 628. The Heterodyne Receiver 
Array CO Line Extragalactic Survey (HERACLES, 
ILerov et al.]l2009| ) C0(2-l) map, obtained at the IRAM 
30-m single-dish telescope, has a beam size of 13.6" and 
rms noise of ^ 22 mK per 2.6 km channel, which 
translates into a la limit on the molecu lar gas surface 
densit y of Y^h^ ^ 3Mq pc~^. Following iSchruba et al.l 
(|2011| ) we assume a constant CO(2-1)/CO(1-0) line ratio 
of 0.7 to estimate I(CO) from the HERACLES data. 

The Berkeley Illinois Maryland Ar ray (BIMA) Survey 
of Nearby Galaxies (BIMA-SONG, iHelfer et al.l [20031) 
CO (1-0) map combines zero spacing single dish data from 
the NRAO 12 m telescope and interferometric BIMA C 
and D array data, resulting in a map with a robust beam 
size of 6.2", and rms noise of 51 mJy beam"-'^ in a 10 km 
s~^ channel or S//^ ^ 13 M© pc~^. 

Finally, the CARMA CO(l-O) map of NGC 628 
() Rahman et al.l i2012) has a spatial resolution of 3.6" and 
rms noise of '^20 mJy beam~^ in a 10 km s~^ chan- 
nel (Sif^ - 5Mq pc-2). We convolve the BIMA-SONG 
and CARMA maps with a Gaussian kernel to match the 
13.6" PSF of the HERACLES map, which corresponds 
to ~ 570 pc at the assumed distance to NGC 0628. The 
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Fig. 1.— Maps of the Ha, from VENGA, and FUV+24^m, from 
Galex and Spitzer-IRAC, SFR surface density (top two), and the 
CO intensity from the HERACLES, BIMA-SONG, and CARMA 
data respectively (bottom three). All maps are convolved to the 
HERACLES beam-size of 13.6" or 570 pc at the distance of NGC 
628 (shown as a white circle in the top right corner of each map). 
Solid lines show the edges of the radial annuli in which we measure 
^CO- Gaps in the maps correspond to regions contaminated by 
foreground MW stars. 



three PSF matched CO maps are shown in Figure [TJ 

In order to account for systematic errors in the ab- 
solute flux calibration of the data-cubes, and in the as- 
sumed C0(2-l) to CO(l-O) conversion factor, we scale 
the maps to match the total CO(l-O) luminosity of the 
HERACLES map in the inner 5 kpc of the galaxy. We 
find scaling factors of 1.36 and 1.36 for the BIMA-SONG 
and CARMA maps respectively. Figure [5] presents both 
the scaled and original flux measurements in each of 
the radial bins over which we will measure Xqo (see 
§3.2) for the three maps. The ~30% difference implied 



by these scaling factors is consistent with the system- 
atic uncertainties in the absolute flux calibration of the 
three datasets, and the observed scatter in the C0(2- 
1 ) to CO (1-0) ra tio across the HERACLES galaxies 
(|Lerov et al.|[2009l ). Assuming that the C0(2-l) to C0(1- 
0) ratio has no significant radial dependance, these cor- 
rections only introduce a scaling in the derived Xco 
factors. Since we are interested in measuring relative 
changes in Xqo within the galaxy, a scaling of this type 
does not affect our results. 

The radial dependance in the C0(2-l) to CO(l-O) ra- 
tio across the disk of spiral galaxies has been studied by 
iLerov et al.l (|2009D by comparing the C0(2-l) emission 
in the HERACLES m aps to the Nobeya ma 45m single- 
di sh CO(l-O) maps of iKuno et all (|2007f ) (see Figure 34 
in ILerov et al.l (j2009() ). Although the scatter seen in the 
line ratio is large, no evidence for a radial dependance 
is seen outside the very central parts of the galaxies 
(i? < 0.05i?25)- Furthermore, as will be discussed in 
§4, the relatively good consistency between the different 
datasets (see also Figure [2|) , and the agreement between 
our measured Xco radial profiles and independent mea- 
surements based on dust modeling indicate that any po- 
tential radial trend in the C0(2-l) to CO(l-O) ratio is at 
a significantly smaller level than the radial trends seen 
in Xco- 

The consistency between the different measurements 
shown in Figure[2|is reassuring. While the BIMA-SONG 
map includes both BIMA interferometric data and 
NOAO 12m single-dish data, the CARMA map is 
constructed without the addition of zero spacing infor- 
mation, and therefore lacks sensitivity on large scales. 
The agreement seen in Figure [2] implies that even in 
the outermost radial bin for which we can measure the 
flux in the CARMA map, the lost large scale extended 
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Fig. 2. — CO intensity as a function of radius for the four radial 
bins adopted in this work. Filled circles correspond to the CO(l- 
0) intensities after applying the scalings described in §2.3. Open 
circles show the raw BIMA-SONG and CARMA measurements 
before scaling them to match the HERACLES absolute flux level. 
Open squares show the raw CO(2-l) intensity in the HERACLES 
map before correcting for the assumed CO(2-l) to CO(l-O) ratio 
of 0.7. Datapoints within each radial bin have been shifted in 
the horizontal direction for clarity. Vertical dotted lines mark the 
edges of each radial bin. The error-bars for each measurement, 
calculated as the r.m.s. per beam in each map divided by the 
square root of the ratio between the area of each radial bin to the 
beam-size, are smaller than the datapoints themselves, and not 
shown. Differences in the profiles are dominated by systematic 
uncertainties at the 20% level. 
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emission does not contribute significantly to the total 
CO surface brightness. 



3. ESTIMATING Xco 

Estimating Xqo (Equation 1) requires an indepen- 
dent measurement of the molecular gas column den- 
sity, or equivalently its surface density {'Eh2), and the 
CO(l-O) line intensity (/(CO)). The latter is directly 
measured from the CO maps described in §2.3. We esti- 
mate T,H2 from the SFR surface density (S.^ ^r) by in- 



verting the molecular gas star formation law (jKennicuttl 
[TOQllKennicutt et al.lf2007l : iBigk 



del et al.l[200l . which we 



parametrize as 



^SFR 



IMgyr-ikpc" 



= A 



^H2 



lOMopc- 



N 



(2) 



This is a more general approach than a ssuming a con- 
stant molecular gas depletion timescale (iSchruba et al.l 
IMllGenzel et ahllMa iMcQuinn et al.ll2012D . which is 
equivalent to the above method for the particular case of 
iV = 1. 

For simplicity, we fix the SFL normalization to 
A = —2.3, which is equivalent to a molecular gas deple- 
tion timescale of 2 Gyr at = 10 pc~^ . This value 
is in good agreement with observations of so lar metallic 



ity spiral galaxies in the local universe (e . g. IBigiel et al 



2008 ; iLero v et al. 20081: iBlanc et al.ll2009t ISchruba et al 



201 Ih . and modifying it introduces a simple scaling in 



the derived Xco values. We report the Xco radial pro 
file for two assumed values for the SFL slope oi N — 1.0 
and N = 1.5, which span the range of plausible slopes al- 
lowed by random and systematic uncertainti es in current 
measurements of the molecular gas SFL (C alzetti et al.l 
120121 : iKennicutt fc Evan£ll2012L and references within) 

It is important to explicitly state the assumptions 
underlying our method to measure Xqo- Mainly, we 
are assuming the existence of a fundamental power- 
law like correlation between 'Sh2 and T.sfr, which 
holds for averaged measurements of these two quan- 
tities over kpc scales. The existence of such cor- 
relation across the disks of spiral g alaxies has ben 
well stablished by previo us studie s ("Kennicu tt et al 



20121 ) 



20071: Bigiel et al. 2008; Bl anc et al.i 2009: Verl ev et al 
201C : rOnodera et al.l2010tlLiu et al.l2011tlRahman et al 



A potential caveat arises from the fact that the ob- 
served SFL in all the above mentioned studies has been 
measured by assuming a constant Xqo factor to trans- 
form CO surface brightness into ^h2 ■ We do not consider 
this to be an important limitation, as reasonable changes 
in Xco (i-e. of the magnitude expected from theoretical 
models for the range of ISM physical properties present 
in galaxies) are not large enough as to break the observed 
correlation betw een TiH2 and Tisfr- On the contrary, a 
recent study bv iNaravanan et al.l ()2012|) shows that us- 
ing a changing Xqq, which depends on the metallicity 
and surface density fo the molecular gas, translates into 
a smaller scatter around the best-fit power-law SFL for 
integrated measurements of galaxies, than when a single 
or a bimodal Xco is used. Therefore, a power-law like 
SFL on kpc scales appears to exist independently of the 
assumption of a constant or changing Xco factor. 



A second assumption made when applying our method 
is that the adopted SFL does not change across the disk 
of the galaxy. In the linear SFL scenario (i.e. when as- 
suming N — 1) this is equivalent to assuming a constant 
depletion timescale across the disk of the galaxy. In the 
linear case, our method cannot break the degeneracy be- 
tween the measured values for Xco and the assumed gas 
depletion timescale. As mentioned above, the absolute 
adopted value only introduces a simple scaling of Xcoj 
which is not important given that we are interested in 
studying relative changes across the disk of the galaxy. 
However, any observed changes in Xco within the galaxy 
could, in principle, be attributed instead to changes in 
the depletion timescale, or equivalently, in the star for- 
mation efficiency (SFE, defined here as the inverse of the 
depletion time). In §4.3 we discuss this possibility, only 
to conclude that it is very unlikely that changes in the de- 
pletion time can explain the observed Xco radial profile, 
and that the observed change in the conversion factor is 
most likely real. In the non-linear case {N = 1.5), the de- 
pletion timescale, or equivalently the SFE, does change 
as a function of the local surface density of molecular 
gas. 

3.1. Measurement of the SFR 

We produce maps of T,sfr using two independent SFR 
indicators. First we use the Ha emission line flux, cor- 
rected for dust extinction using the Balmer decrement. 
The methods used to measure dust corrected Ha fluxes 
and to correct the emission line map for the contribution 
from diffuse ionized gas ( DIG) are analog to the ones 
described in IBlanc et all ()2009[ ). An updated descrip- 
tion of the methods used, including slight differences in 
the methodology used for flux calibration and construc- 
tion of the IFU data-cubes will be described in an up- 
coming public a tion. We use the Ha SFR c alibration of 
iMurphv et all (120111) andlHao et al.l (120 111) . taken from 
the com pilation in Kennicutt fc Eyana ( 2012l ). which as- 
sumes a iKroupa fc Weidneii (12003^ initial- mass- function 
(IMF). Note that this d i ffers f rom the ISalpeteJ ()1955[ ) 
IMF used in IBlanc et all (pOOOt ). 

The second SFR indicator used is a linear combina- 
ti on of FUV a nd 24 ^m flux. We us e the calibration 
of ILerov et all (|2012D . which assumes a IChabrieri (l200l 
IMF. Calibrations based on the Kr oupa and Chabrier 
IMFs yield nearly identical SFRs (jChomiuk fc PovichI 
1201 If ), which are typically ^-^30% lower than those ob- 
tained assuming a Salpeter IMF. Since Ha and FUV 
photons trace star formation over different timescales 
(roughly 10 and 100 Myr respectively), systematic differ- 
ences between the two independent estimates might arise 
as a consequence of differences in star formation history. 
Analogously to what was done with the CO maps, we 
attempt to remove any systematic differences in the flux 
calibration of the different datasets and the SFR calibra- 
tions used, by matching the total SFR within the central 
5 kpc of the galaxy. This translates in us multiplying 
the 'SsFR,FUV+24:fj.m map by a factor of 0.8 to match the 
'^SFR,Ha map. Again, this correction does not affect the 
relative differences in Xco which we are trying to mea- 
sure, but only introduces a scaling of the derived values 
across the whole galaxy. Both maps are convolved with 
a Gaussian kernel in order to match the common 13.6" 
FWHM PSF, and are presented in Figure [TJ 
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3.2. Calculation of Xqq 

We use the PSF matched SFR maps to measure 
^SFR,FUV+2ifj,m, ^SFR,Ha, and I (CO) in ehiptical an- 
nuH of constant galactocentric radius. We chose an an- 
nuli width of 40" which roughly corresponds to three 
times the spatial FWHM of the maps, and a physical 
scale of 1.68 kpc at the adopted distance to NGC 628 
(see Figure [T]). The sensitivity of the BIMA-SONG and 
CARMA maps only allows the measurement of I{CO) 
out to 5 kpc (~ 120") which corresponds to the in- 
ner three annuli, while the HERACLES map is able 
to provide a reliable measurement in a fourth bin, out 
to a galactocentric d istance of 7 kpc (~ 170"). While 
ISchruba et all (poTTI) was able to measure CO emission 
using the HERACLES map of NGC 628 out to larger 
radii (~10 kpc) using a stacking technique, here we are 
also limited by the coverage of the VENGA IFU obser- 
vations which only reach out to ~7.5 kpc. We use the 
extended stacking measurements of lSchruba et al.l (|2011[ ) 
separately in the next section to extend our measure- 
ments to larger radii. Having measured these three quan- 
tities {'£.SFR,FUV+2iiira, '^SFR,Ha, and /(CO)), we USC 

Equation 2 to obtain E^f2, and we transform it to units 
of column density. Finally we calculate Xqq for each 
radial bin by applying Equation 1. The central radius 
and measured Xqq values for each annuli are reported 
in Table m 

The uncertainty in the measured Xqq values is domi- 
nated by the intrinsic scatter in the molecular SFL. A se- 
ries of studies have shown that this quantity depends on 
the physical scales ov er which the observed surface den- 
sities arc integrated ("Verle v et al.l I2010l: lOnodera et al.l 
IMot rSchruba ct al. 2011; Liu et al.ll20'Tll ). For the large 
physical scales considered in this work ('^ 4 kpc given 
the mean area of our radial bins), we assume a value 
of ~ 0.15 dex for the scatter, which we propagate into 
the uncertainty in Xqo- Thanks to the high S/N of the 
datasets and the large areas over which we are integrat- 
ing, the contribution from random errors in Yisfr, which 
includes photometric errors in the Ha and Hj3 (the lat- 
ter entering in the B aimer decrement dust extinction cor- 
rection), and in the FUV and 24 /im fluxes, as well as 
from noise in the CO maps, is negligible (a factor of ~ 5 
smaller) compared to the scatter in the molecular SFL. 

4. THE Xco RADIAL PROFILE 

4.1. Consistency Between Different Datasets 

In Figure |3] we present the measured Xqo factor as 
a function of galactocentric radius for the six combina- 
tions of datasets (two SFR indicators and three CO maps 
from different telescopes) , and two different assumptions 
for the molecular SFL slope (A^ = 1.0 and 1.5). At aU 
radii, all combinations of datasets are consistent with 
each other within the la formal error bars, and show the 
same systematic trend of increasing Xco with radius. 
This implies that our results are robust against potential 
systematics associated with the adopted SFR indicator, 
and the nature of the CO maps utilized (single-dish vs. 
interferometer, CO(l-O) vs. C0(2-l)). The consistency 
between different datasets also confirms that photometric 
errors are not a significant source of uncertainty, when 
compared to the systematic errors associated with the 
scatter in the SFL. 



As can be seen in Figure (TJ the agreement between the 
Ha and FUV-f 24/xm SFR maps is good. After perform- 
ing the scaling described in §3.1, the ratio of the two dif- 
ferent SFRs measured in the inner four radial bins shows 
a r.m.s. scatter of 13%. When analyzing the morphol- 
ogy of these two SFR tracers we notice a strong similar- 
ity. All individual star forming regions picked up by one 
method are also recovered by the other. This is impor- 
tant as it implies that the dust corrected Ha method is 
not missing a significant fraction of completely obscured 
star formation which should be recovered by dust emis- 
sion at 24/im. 

While both maps trace the same structures, a more de- 
tailed inspection shows that the FUV-|-24y^m SFR map 
does seems to have a more diffuse morphology than the 
Ha SFR map. In particular more diffuse FUV-|-24/^m 
emission is seen in the inter- arm regions of the galaxy, 
while the Ha emission seems to be more concentrated 
in the arms. This is consistent with the typical star 
formation timescales traced by these two methods. Ha 
traces young stars (< 10 Myr), which are seen close to 
their birthplaces, while both FUV and 24/im emission 
receive a significant contribution from older stars (< 100 
Myr) which are more homogeneously distributed across 
the galactic disk, and not necessarily associated with the 
spiral arms. 

The morphological agreement between the different 
CO maps is also very good. After applying the scal- 
ings described in §2.3 the CO luminosities in the inner 
3 radial bins show a typical r.m.s. scatter of 20% (see 
Figure [5]). The agreement justifies our adoption of a 
single CO (2-l)/(l-0) line ratio across the galaxy. It is 
interesting to notice that while that both the CO emis- 
sion and the SFR trace similar morphological structures 
across the galaxy, that is mainly the grand design spi- 
ral arms of NGC 628, a detailed comparison shows that 
there is not a one to one correlation between individual 
CO bright regions and regions which are bright in Ha 
or FUV-|-24//m. This refiects the temporal dependance 
of the star formation process, in which recently formed 
star clusters disrupt and dissociate their molecular gas 
birth clouds by injecting both radiative and mechanical 
energy into their surrounding ISM. As this process takes 
place, and individual star forming region evolves from a 
CO bright molecular complex into a 24/im bright, highly 
obscured, star forming region, which eventually dissoci- 
ates and removes the bulk of the molecular gas and dust 
surrounding it, becoming a Ha and FUV bright region 
in the process. This temporal evolution translates into 
an spatial offset between SFR and molecular gas tracers, 
and is the main cause behind the need to average over 
large kpc scales in order to do a measurement like the 
one p resented in this work (Schr uba et al .''2010': 'Mador^ 
120101: iFeldmann et al.l [201 2b : .Kennicutt fc Evans 2012,1 

In order quantify the previous statements regarding the 
robustness of our measurement we calculate the standard 
deviation of all Xqq measurements in each radial bin, 
and find it to range from a few percent to a ~ 25% of the 
measured values in the worst case (outermost bin). The 
mean r.m.s. across all four bins corresponds to ~ 15% 
of the measured values. This dispersion is caused by a 
combination of different factors including systematic er- 
rors in the adopted CO(2-1)/CO(1-0) ratio, a possible 
systematic lack of sensitivity to low surface brightness 
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Fig. 3. — Radial profile of Xco in NGC 628. The top and bottom panels present results for a SFL slope of N = 1.0 and A'^ = 1.5 
respectively. Measurements using the HERACLES, BIMA-SONG, and CARMA maps are shown in green, blue, and red respectively. 
Measurements based on Ha and FUV+24/im SFRs are shown as large color squares and circles respectively. Datapoints within each radial 
bin have been shifted in the horizontal direction for clarity. Yellow stars show the average of these measurements in each radial bin. Vertical 
dotted lines mark the edges of each radial bin, and the horizontal dotted line marks the MW Xco factor of Pineda ct al. (2010). Also 
shown are measurements from Spitzer+ Herschel dust SE P modelling (Karin S andstrom private communication, small open circles), and 
the results of applying our method to the stacked data in lSchruba et aLl II2011I) (small filled circles). 
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emission in the interferometric CO maps, and potential 
variations in the ratio between the Ha and FUV+24/L(m 
SFRs which might arise as the consequence of changes in 
the stellar populations and the properties of dust across 
the disk of the galaxy. The magnitude of these systematic 
deviations is smaller than the formal 0.15 dex 35%) 
uncertainty adopted. Nevertheless, we add this disper- 
sion in quadrature to our formal errors in order to take 
into account the effects of theses systematics in our error 
budget. 

4.2. Dependance of Xqq on the Assumed SFL Slope. 

The dotted horizontal line in Figure[3]marks the canon- 
ical Milky Way Xco fact o r of 2 .3 x lO^o cwr^{¥. km 
s^^)~^ from I Pineda et al.l (120 lOD. which is ado pted as 
the preferred value in lKennicutt fc Evans! ()2012l ). While 
at radii > 2 kpc the Xco factor in NGC 628 is consis- 
tent with the MW value, there is an evident drop in the 
central 2 kpc of the galaxy, where the measured Xqq is 
a factor of 2 lower than the MW value. 

The drop in the radial profile of Xqq towards the cen- 
tral regions is observed independently of the assumed 
SFL slope. A non-linear SFL with TV > 1 implies lower 
H2 surface densities, and therefore lower derived Xco 
values for high SFR surface density regions than a lin- 
ear relation. Since the SFR in NGC 628, and in most 
disk galaxies, decreases with radius, a non-linear slope 
implies an even larger drop in Xqq towards the central 
regions. In any case, the magnitude of the changes in 
Xco between the = 1.0 and 1.5 cases is smaller than 
the changes in Xqq observed as a function of radius 
across the galaxy. This implies that potential changes 
in the SFL slope within the galaxy would not affect our 
results significantly. While a sub-linear (A^ < 1) slope 
could reduce the magnitude of the observed drop in Xco 
towards the central regions, this scenario is not cur- 
rently supported by obs ervations (jCalzetti et al.l 120121 : 
iKennicutt fc Evansll2012l and references within). 

This dependance with the SFL slope is important, as 
it implies that the observed radial trend of Xco is not 
a consequence of changes in the star formation efficiency 
with radius. Adopting a super-linear slope for the SFL is 
effectively testing the case in which the efficiency of star 
formation is enhanced in higher density regions, and as 
stated above, this effect only boosts the decrease of Xco 
at small radii, although not significantly. We further 
discuss this subject in the next section. 

4.3. Xco Gradient and Comparison to Other 
Measurements 

Also shown as open circles in Figure [3] are the mea- 
surements of Xco in NGC 628 from dust SED modelling 
(Karin Sandstrom private communication). We only plot 
their data-points with an uncer tainty < 0.4 dex in Xco- 
This study uses the method of iLerov et al.l ()20lH ) to si- 
multaneously model the dust mass surface density and 
the dust-to-gas ratio across the disks of a sample of 
nearby spiral galaxies with far-I R SED measurements 
from the SINGS and KINGFISH (jKennicutt et al.|[20Tll) 
surveys. The dust mass surface density, in combination 
with the dust-to-gas ratio, provides an estimate of the 
total gas surface density. Subtraction of the atomic com- 
ponent (measured from HI 21 cm maps) yields the surface 



density of the molecular component. This technique is 
completely independent of the method adopted in this 
work. The data-points from the dust SED modelling 
have a typical uncertainty of 0.2 dex at i? < 7 kpc, hence 
they agree with our measurements across the whole range 
in radii sampled by our data, and significantly show the 
same trend of decreasing Xco towards the central re- 
gions of NGC 628. No scaling has been applied to the 
measurements in Figure [Sj so the agreement is not only 
in the shape of the Xco radial profile but also in its 
absolute value. 

A s imilar analysis was conducted by lAniano et al.l 
(|20ll) for NGC 628 and NGC 6946. While a drop in 
Xco towards the center of the galaxy is evidently seen in 
their data for NGC 6946, the authors do not claim the de- 
tection of a changing Xco profile in NGC 628. Inspect- 
ing their Figure 4, its is evident that the results of their 
measurements of Xco in NGC 628 are subject to strong 
systematic uncertainties, and show large changes (from 
a decreasing to a flat, or even decreasing radial trend) 
depending on which subset of the Spitzer+Herschel data 
is used and the spatial resolution at which the modeling 
is conducted. This is not the case for NGC 6946 where 
different datasets at different resolutions yield consistent 
results. Given the large systematic uncertainties in their 
measurement of Xco for NGC 628, we refrein from con- 
ducting a detailed comparison to their results. 

In §3, we discussed the inherit degeneracy in our 
method between the SFL normalization constant (or 
equivalently the gas depletion time when N = I), and 
the measured values of Xco- A basic assumption in this 
work is that of a constant SFL across the disk of NGC 
628. In order to explain the data in Figure [3] solely as a 
variation in the normalization of the SFL would require 
the depletion timescale to be longer, or the star forma- 
tion efficiency to be lower, by a factor of 2 in the central 
(denser) regions than in the outer disk. This is not in line 
with expectations from measurements of the SEE across 
different environments in galaxies in the local universe. A 
series of studies have shown that luminous infrared galax- 
ies (LIRGS), ultra-luminous infrared galaxies (ULIRGS), 
and sub-mm galaxies (SMGs), in which the molecular gas 
surface density is one to two orders of magnitude larger 
than in the disks of normal spirals, show typical deple- 
tion timescales which are factors of 3-4 shorter than nor- 
mal spiral galaxies. The difference becomes even larger 
(factors of 4-10 shorter) if differences in the Xc o fac- 
tor in these systems are taken into accoun t (Daddi et al.l 
[2010t IGenzefeFalllMl IGarci'a-Burilioet al...2012l) . Al- 
though the central region (< 2 kpc) of NGC 628 is only 
factors of a few denser than the outer disk, it would be 
surprising to find a reversal from the global trend seen 
in the SFE towards denser environments. Furthermore, 
the consistency between the results from this work and 
the independent dust SED modelling method, which is 
not subject to such a degeneracy between Xco and the 
depletion timescale, is very encouraging, and gives us 
confidence that our working assumption is valid within 
the range of current uncertainties. 

If we model the Xco profile in the inner 7 kpc of the 
galaxy as a simple gradient of the form: 

log(Xco) = log(Xco.o) + Alog(Xco) x R (3) 



Xco Gradient in NGC 628 



9 



TABLE 1 
NGC 628 Xco Radial Profile 





N 




Xco - Ha^' 






Xco - FUV+24/xm=' 










HERACLES 


BIMA 


CARMA 


HERACLES 


BIMA 


CARMA 




0.84 


1.0 


1.4±0.5 


1.0±0.4 


1.3±0.5 


1.6±0.6 


1.2±0.4 


1.5±0.5 


1.3±0.5 


2.52 




1.7±0.6 


1.9±0.7 


1.8±0.6 


1.7±0.6 


1.9±0.7 


1.8±0.6 


1.8±0.6 


4.20 




2.1±0.8 


2.6±0.9 


2.1±0.8 


1.9±0.7 


2.4±0.8 


1.9±0.7 


2.2±0.8 


5.88 




1.7±0.8 






2.5±1.0 






2.1±0.9 


0.84 


1.5 


1.3±0.5 


1.0±0.4 


1.2±0.4 


1.4±0.5 


1.0±0.4 


1.3±0.5 


1.2±0.4 


2.52 




1.7±0.6 


1.8±0.6 


1.8±0.6 


1.6±0.6 


1.7±0.6 


1.7±0.6 


1.7±0.6 


4.20 




2.0±0.7 


2.5±0.9 


2.0±0.7 


1.7±0.6 


2.2±0.8 


1.8±0.6 


2.0±0.8 


5.88 




2.6±1.0 






3.0±1.1 






2.8±1.0 



In units of 10^" cm~^(K km s"^)"^ 

a fit to the average of the six measurements in each 
radial bin (two measurements in the outer bin), which 
are shown as stars in Figure [3l yields Xco gradients of 
Alog(Xco) = 0.06 ± 0.02 dex kpc"' for N = 1, and 
Alog(Xco) = 0.08 ± 0.02 dex kpc'^ for N = 1.5. 

The best-fit central Xco values are Xco.o ~ 1.3±0.2 x 
10^0 cm-i(K km a'^)'^ for iV = 1 and Xco,o = 1-1 ± 
0.2 X lO^o cm-i(K km s-^)-^ for the N = 1.5 case. Our 
data imply that Xco in the galactic center of NGC 628 
is about a factor of two higher than the typical values of 
Xco,o = 0.1 - 0.5 X 10^0 cm-i(K km s measured 
in the MW Galac tic Center (e.g. fSodroski et all 119951: 
IStrong et a l. 2004; Oka et al.|[l998[) . 

Inspecting Figure [31 particularly for the A'^ = 1 case, it 
is not evident if the radial profile of Xco follows a linear 
gradient, or if it is fairly flat across the disk of the galaxy, 
and only drops in the central regions. Our measurements 
are limited at i? > 7 kpc by the sensitivity of the HER- 
ACLES CO map, and the coverage of the VENGA IFU 
data-cube. The dust SED modelling datapoints go out 
to ^ 10 kpc, but they show a very large scatter, and at 
such large radii, the uncertainty in their measurements 
increases sharply, so it is not easy to draw conclusions 
from these data-points regarding the behavior of Xco at 
large radii. 

In an attempt to extend our measurements into the 
outer disk of the galaxy, we use our method to estimate 
Xco on the stacked CO and FUV-|-24/xm radia l profiles 
of NGC 628 presented in iSchruba et al.l (|20lH) . In this 
work, the authors use the 21 cm HI velocity field to pre- 
dict the exact wavelength of the CO line at every posi- 
tion in the HERACLES data-cube. This information is 
then used to register and stack the spectra of all spa- 
tial resolution elements in 15" w ide radial bins across 
the galaxy. Using this technique, ISchruba et al.l (|2011[ ) 
is able to significantly measure liCO) out to ^ 10 kpc. 
They also measure the FUV+24/im SFR profile using the 
same maps used in this work. We show these measure- 
ments as filled black circles in Figure [31 Error-bars are 
not shown for clarity, but they are also dominated by the 
intrinsic scatter in the SFL, and are of the same order 
of magnitude as the error-bars for the measurements in 
the inner regions of the galaxy (color filled circles and 
squares). The results of this exercise are inconclusive. 
While in the = 1 case the value of Xco seems to 
flatten around the canonical MW value, if we assume 
N = 1.5 Xco such hints for a flattening in the outer 



regions are less obvious, but still present. 

As will be further discussed in §5, NGC 628 shows a 
metallicity gradient desc ribed by a fairly constant s lope 
out to -1 i?25 (--13 kpc 'Rosale s-Ortega et aT]|201lD . If 
Xco scales as some power of the metal abundance {Z), 
we would expect it to keep on increasing out to a similar 
radius. On the other hand there are reasons why a flat- 
tening in the radial profile might be expected. For exam- 
ple, if Xco increases with a de creasing molecu lar cloud 
surface density, as proposed bv lNaravanan et~a l. (2012), 
and the GMCs in the outer disk have uniform proper- 
ties but become denser towards the central regions, a 
certain level of flattening in the Xco profile might oc- 
cur, although the metallicity dependance discussed above 
should still be present. We will explore the role of the 
molecular gas smface density and the UV radiation field 
in §6, and will discuss this possibility in more detail. 

For all the results presented in the following sections, 
we adopt the iV = 1.0 case. We have checked that 
changing the slope to = 1.5 does not significantly 
affect the rest of the results presented below. 

5. RELATING THE Xco AND METALLICITY GRADIENTS 

As discussed in §1, both observations of nearby and 
high redshift galaxie s (lArimoto et al.1 119961: [Israeli il997l: 
Bolatto et all [20081: [Lerov et al.[ [20111: [Schruba et al.l 
20121: [Genzel et al.[ [2012[) as well as theoretical 
models ([Krumholz et al. l |2011[ : [Naravanan eFall [2011: 
[Feldmann et al.|[2012a[ ) iiidicate that Xco is a decreasing 
function of metallicity. The main physical process behind 
this dependance is the role that dust extinction, which is 
directly linked with metallicity, plays in the dissociation 
of CO molecules. In photo-dissociation regions on the 
edges of molecular clouds, as the metallicity and with 
it the dust extinction decreases, the CO to C+ transi- 
tion layer moves inwards, leaving behind large envelopes 
of "CO dark" gas where hydr ogen is still in molecu- 
lar form thanks to self- shielding ( Malonev fc Blacl3[1988[ : 
[Bolatto et al[[l999l: [Glover et al.|[2oTo[ r This translates 
into higher CO to H2 conversion factors for lower metalli- 
cally environments. While at the relatively high metallic- 
ities found across the disk of NGC 628 we do not expect 
to see extreme cases of completely "CO dark" molecular 
clouds (which might be present in extreme lo w metallic- 
ity galaxies like I Zw 18. [Schruba et al.|[2012[ ). variations 
in the CO brightness of the diffuse lower density (and 
lower extinction) envelopes of molecular clouds as func- 
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tion of metallicity are possible ()Feldniann et al.ll2012af) . 

Since most disk galaxies show radial metallicity gradi- 
ents, wi th higher metal abunda nces towards the central 
regions (|Moustakas et al.]l2010L and references therein), 
it would be natural to attempt to relate metallicity gra- 
dients with gradients in Xqq- In this section we study 
the relation between the observed radial profile of Xco in 
NGC 628, and its metallicity distribution. We measure 
the metallicity gradient using strong nebular emission 
lines in the VENGA IFU data, and study how the rela- 
tion between Xqq f^nd metallicity across the disk of the 
galaxy compares to that seen across galaxies in the local 
universe, and to the predictions of theoretical models. 

5.1. The Metallicity Gradient in NGC 628 

The VENGA IFU data-cube provides maps of 
many strong nebular emission lines typically used for 
metallicity diagn o stics. As discussed in d e tail i n 
iKewlev fc EllisonI (|2008[ ) and iMoustakas et all (|2010[ ). 
different strong line methods used to measure metallic- 
ity can show large discrepancies in the derived abun- 
dances. Particularly, methods calibrated against theoret- 
ical photo-ionization models typically yield higher metal- 
licities than methods calibrated against samples of HII 
regions with direct electron temperature measurements 
(a.k.a. empirically calibrated methods). Differences as 
large as 0.6 dex are seen between these two families 
of methods. On the other hand, relative differences in 
metallicity from one system to another do not change 
drastically if different methods are applied, allowing 
meaningful conclusions to be drawn from comparing dif- 
ferent samples, as long as a single method is used, or mea- 
surements using different methods a re properly trans- 
forme d to a common metallicity scale (|Kewlev fc EllisonI 
120081 ). In order to study the metallicity dependance of 
Xco across NGC 628, and compare it to previous mea- 
surements in the literature, we must not only choose a re- 
liable strong-line abundance indicator, but we must also 
put ours and the literature measurements on a common 
metallicity scale. 

Our preferr ed metallicity calibration is the 7V202 
method of .Kewlev fc Dopital (|2002j ). which estimates the 
oxygen abundance from the [NI]A6584/[011]A3727 ratio, 
and is calibrated against photo-ionization models. We 
have a series of reasons to favor this calibration. Mainly, 
the N202 indicator is highly independent of the ioniza- 
tion parameter, is single valued at all metallicities, and 
the secondary production of nitrogen makes this indica- 
tor very sensitive to metallicity changes at high metallici- 
ties like the ones we expect i n massive spiral galaxies like 
NGC 628, (see discussion in IKewlev fc Dopital [200l . 

On the other hand, the metallicity compilations for 
the systems \ yith m eas ured Xgo values reported in 
iBolatto et all (|2008f ) and iLerov et al.l (j2011[ ) are largely 
based on direct electron temperature measurem ents. The 
metallicities reported in iSchruba et al.l (j2012|) are on an 
intermediate scale, since the authors averaged the val- 
ues from the two indicators used by IMoustakas et al.l 
(pblO) : t he photo-ionization model cali brated R23 metal- 
licities of iKobulnickv fc KewlevI (|2004[ ) (hereafter R23 - 
KK04) and the empiric ally calibrated R23 method of 
iPilvugin fc ThuanI (|2005h (hereafter R23 - PT05). 

We decide to apply a simple offset in metallicity to 
both our N202 measured values in NGC 628 and the re- 




FlG. 4.— Spatial N202 metallicity distribution in NGC 628 mea- 
sured from the VENGA data-cube, offset to match the "empirical 
method" metaUicity scale. The top panel shows a two-dimensional 
map of 12-f log{0/H). Thick contours marks regions above the Ha 
surface brightness cut where the nebular emission is dominated by 
HII regions. The bottom panel shows the radial distribution in 
metallicity of all spaxels inside these regions (black dots). Red cir- 
cles and error-bars mark the median and standard deviation of all 
spaxels in 0.5 kpc radial bins. Black error-bars show the median 
measurement error in each radial bin, and green error-bars show 
the 0.1 dex intrinsic scatter associated with the N202 calibration 



ported metaUicity values in ISchruba et al.l (j2012l ). to put 
them on a commo n scale with the dire ct me thod metal- 
licities reported in iBolatto et al.l (|2008f) and ILerov et al.l 
((20T1I) . To do so, we measured the metallicity in the 
NGC 628 data-cube using the three methods mentioned 
above (iV202, R23 - KKOA, and R23 - PTOS), and mea- 
sure the mean offset in the recovered metallicity val- 
ues between different methods for all spaxels above a 
Ha surface brightness cut chosen to trace HII regions 
Fl . As expected, both the N202 and the average be- 
tween the R23 - K K04, and R23 - PTOb methods (i.e. 
the method used in ISchruba et al.l (|2012l )) yield higher 
metallicities than R23 — PT05 alone, which is calibrated 
against direct method measurements on HII regions and 
should fall in a sim ilar s cale to the values reported in 
IBolatto et all (f2008l) and ILerov et al.l (|20ni ). We find 
and apply me an offsets of 0.66 dex and 0.35 dex to the 
N202 and Sc hruba et al.l (j2012[ ) metallicities before do- 
ing any comparisons. These values are in good agreement 
with the 0.6 de x offset found between R 23 — KKOA, and 
R23 - PT05 bv IMoustakas et al.l (|2010l ). 

Note that we decided to apply a single offset in 
metallicity instead of u sing the conversion relations of 
IKewlev fc EllisonI (|2008l ). This is done deliberately, as 

^■^ All these strong-line methods have been calibrated either 
empirically against direct-method metallicities of HII regions, or 
photo-ionization models with densities and ionization parameters 
typical of HII regions, therefore they cannot be applied to low sur- 
face brightness regions in which the nebular emission is dominated 
by the diffuse ionized gas. 
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the non-linear nature of these conversions imply that the 
measured metallicity gradient would change its value de- 
pending on the adopted scale. Also, these conversions 
have been calibrated against integrated measurements of 
galaxies, and it is not clear how valid it would be to apply 
them to spatially resolved regions within a single object. 
This issue will be the subject of a future publication in 
which we will analyze the impact of different methods 
on the observed metallicity distributions in the VENGA 
galaxies. For the purpose of this work, we prefer to con- 
serve the shape of the gradient as measured by the N202 
method which we consider more robust than other formu- 
lations. In any case, using the lKewlev fc EllisonI ([2008h 
conversions, or using R23 — PT05 directly on our data in- 
stead of N202 does not significantly change our results. 

Figure [4] presents the spatial distribution in metallic- 
ity of NGC 628. The top panel shows a two-dimensional 
map of 12-|-log(0/H), where thick contours marks re- 
gions above the Ha surface brightness cut where the neb- 
ular emission is dominated by HII regions. The bottom 
panel shows the distribution of metallicity as a function 
of galactocentric radius for data-cube spaxels inside these 
regions. Black error-bars show the median measurement 
errors in each radial bin, red error-bars mark the stan- 
dard deviation for all spaxels in each bin, and green error- 
bars show the 0.1 de x intrinsic scatter asso ciated with 
the iV202 calibration fKewl ev fc Ellison|[2C)080 . We mea- 
sure a gradient of Alog{OH) = 0.036 ± 0.002 dex kpc'^ 
in good agreement with the IFU measurements of 
iRosales-Ortega et all (|2011[ ). 

5.2. ^co ^s- Metallicity 

In Figure [SI red stars show the Xqo values measured 
in each radial annuli in NGC 628, as a function of the 
metallicity at the radius corresponding to the center of 
each annuli (evaluated from the metallicity gradient mea- 
sured in §5.1). Also shown are the data of Bolatto et al. 
((2008I) . iLerov et all (l20Tll) . and Schruba et al. (20120 
as well as predi c tions from the theoretical models of 
iNaravanan et all (j2012| ) and lFeldmann et all (j2012al ). 

Across the disk of NGC 628, Xco increases towards 
lower metallicity regions. This is expected from the 
positive and negative gradients measured for Xco and 
the oxygen abundance respectively (§4 and §5.1). This 
behavior is consistent with integrate d measurem e nts o f 
galaxies in the local universe. Both iLerov et al.l (poTTI ) 
and lSchruba et al.l(|2012| ) see an increase in Xco towards 
low metallicities by conducting integrated measurements 
across galaxies using dust mass modeling in the first case, 
and a similar method to the one used here in the later. 
The Xqo values in ferred from CMC virial masses in 
iBolatto et al.l (l2008f ) do not show an obvious trend with 
metallicity, but are consistent with other studies given 
the large scatter in t he observed re lation. Also, as dis- 
cussed in the o riginal IBolatto et al.l ([2008) paper, and in 
iSchruba et al.l (2012), estimates of Xco using virial mass 
measurements from CO observations most likely miss the 
"CO dark" envelopes of molecular clouds in low metallic- 
ity galaxies, as they can only trace mass enclosed in the 
CO emitting region. This could translate in an underes- 

The data-points of lSchruba et al] 120121) have been scaled up 
by a factor of 1.11 to account for the difference in the assumed 
depletion timescale (2.0 vs 1.8 Gyr) 



timati on of Xco at low metallicity in the IBolatto et al.l 
(|2008l) data. The same effect could expla in the constancy 
of Xco seen bv lRosolowskv et al.l ()2003D across 45 GMCs 
in M33, even in the presence of a 0.8 dex change in metal- 
licity within their sample. 

The NGC 628 datapoints are consistent with previ- 
ous observations, specially considering the large scatter 
seen in the Xco metallicity relation (~ 0.5 dex). The 
limited dynamic range in metallicity that we can sam- 
ple across the inner 7 kpc of a single galaxy prevents us 
from being able to put strong constraints in the slope of 
the Xco-metallicity relation. Mod eling this relati o n as a 
power- law of the form Xco k Z", ISchruba et all (|20l2l ) 
measures a slope of a = —2.0 ± 0.4 for the full HER- 
ACLES sample after rejectin g starb urst galaxies (green 
solid line in Figure |5|), while llsraell fl997) finds a slope 
of a = —2.7 ± 0.3 in their study of dust emission in local 
dwarfs. A linear fit to our data points yields a slope 
of a = —1.1 ±1.6, shallower than previous mea sure- 
ments but stil l cons istent with both llsraell ()1997[ ) and 
ISchruba et all (I201I) . A gain, given the limited dynamic 
range in metallicity of our data, and the large scatter in 
the Xco-nietallicity relation, detailed comparisons are 
difficult. 

The slope of the observed relation in NGC 628 is also 
in agreement with expectations from theoretical models. 
At fixed molecular g as su rface density, the simulations 
in INaravanan et al.l ()2012| ) predict a linear dependance 
of Xco w ith metallicity ( i. e. a = -1.00, and the 
study by iFeldmann et al.l (|2012aD predicts slopes in 
the —0.5 to —0.8 range, depending on the treatment 
given to the sub-grid gas dynamics in their simulations. 
Regarding the normalization of the observed relation, 
our data shows a rema rkably good agreement with the 
INaravanan et al.l (|2012f ) model for an assumed molecular 
surface d ensity of Sffa = 100 M mpc~^. On the other 
hand, the IFeldmann et al.l (|2012aD model,a = -2.7 ±0.3 
in which Xco is independent of the molecular gas 
surface density, predicts values that are a factor of 4 
higher than our observations. In the following section 
we discuss the characteristics of these models in more 
detail, paying particular attention to the discrepancies 
which arise from the different treatments given to the 
role of the local gas surface density at setting the Xco 
conversion factor. 



6. THE ROLE OF GAS DENSITY AND THE UV 
RADIATION FIELD AT SETTING Xco 

In the previous section we explored the role of metal- 
licity in setting the value of Xco , and discussed how well 
observations are described by the dependance on metal- 
licity predicted by current theoretical models of CO emis- 
sion in molecular gas. Here, we explore the dependance 
of Xco on other physical properties like the surface den- 
sity of molecular gas, and the local UV radiation field, 
which at the same time, affect the temperature and de- 
gree of turbulence in the ISM. Informed by the predic- 
tions of theoretical models, we examine these potential 

Note that in their model this exponent changes if one 
asumes constant CO brightness instead of constant because 
of the metallicity depen dance of Xco (see Equations 6 and 8 in 
INaravanan et al] I|2012I 1'I 
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Fig. 5. — ^CO factor as a functi on of metallic i ty (12+1ok(0/H)). Red st ars c orrespond to the four radial bins in NGC 628 analyzed 
in this work. The measurements of lBolatto et al.1 II2008 I). ILeroy et al.l II2 011I ). and ISchruba et al.l ||2012|) are shown as orange, purple, and 
green circles. The solid and da shed green lin es show the best-fit relation to the clean HERACLES sample in Schruba ct al. (2012). The 
theoretical models of FNaravanan et al.l 120121) and H'teldm ann ct al. ( 2012a|) are also shown. For the Narayanan ct a l. (2012) model, the 
solid blue line show the expected scaling of Xqq with metallicity for a constant J^h2 = 100 Mqpc~^, and the blue dashed lines mark the 
scatter in Xqq seen across their simulated galaxies. Blue stars mark the p redictions from the model taking into account the differences 
in gas density for each radial bin in NGC 628. For the'Feldma nn et al.l ((2012a) model the two solid cyan lines mark, from top to bottom, 
their "constant line-width" and "virialized" models at 4 kpc resolution. Dashed lines mark factors of 2 above and below the two models 
respectively, which the authors state bound the expected scatter in Xqq for their models. 



trends in the NGC 628 data, and explore if other param- 
eters, beyond the metallicity, are important at setting 
the observed Xqo radial profile. 

6.1. Further Comparison to Theory: Gas Density 

In the theoretical model of iFeldmann et aTl ()2012aD . 
Xco is a strong function of the molecular gas surface 
density of a given GMC, with the conversion factor 
growing away from its minimum towards both low and 
high densities. This behavior is a consequence of two 
different effects taking place in the low and high den- 
sity regimes. At high column densities the numerator 
in Equation 1 increases linearly, while the denominator 
(i.e. the CO intensity) is either constant or increases 
sub-linearly (as in the virialized cloud approximation for 

which I (CO) (X Av (X y/N(ll^). This effect is sup- 
ported by spati ally resolved observati ons of GMCs in 
the Milky Way (jHeiderman et al.ll2010[ ). Therefore Xqo 
increases towards high column densities. On the other 
hand, towards lower N{H2), the CO abundance drops 
due to CO dissociation as the dust extinction (which is 
proportional to the column density) decreases, giving rise 
to large "CO dark" envelopes and also implying an in- 
crease of ^co- Both the location of the minimum, and 
the value of Xco at the minimum are strong functions 



of the metallicity in the model. 

When attempting to study these trends on kpc scales 
within or across galaxies, we must take into account the 
fact that we measure averages for ensembles of molecular 
clo uds with a distribution in their physical properties. In 
the IFeldmann et al.l (|2012aD model, this average washes 
away the dependance on column density seen on small 
scales, and the authors predict almost no variation in 
Xco as a function of II2 surface density. It is important 
to point out that in this model the molecular gas tem- 
perature is fixed to 10 K in the calculation of /(CO), and 
the authors either impose a constant CO line-width, or 
assume virialized clouds (i.e. Av cx ^JN{H2)). 

As discussed in the p reviou s section, the m odels of 
iNaravanan et al.l ()2012[ ) and IFeldmann et al.l (j2012al ) 
predict similar behaviors for Xco as a function of metal- 
licity. On the other han d, the simulations conducted by 
INaravanan et al.l ()2012D do not assume a fixed tempera- 
ture for the molecular component of the ISM, and instead 
include relevant heating and cooling calculations, which 
allow the authors to trace the molecular gas temperature. 
Furthermore, for molecular clouds which are larger than 
the spatial resolution of the simulation, they use the sim- 
ulations themselves to follow the dynamical state of the 
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molecular gas (i.e. the gas velocity dispersion sets the 
CO line- width). For sub- resolution clouds they assume 
virialization. 

The ability to follow the temperature and d ynamical 
state of the gas, allows ;Naravanan et al.l (|2012 ^ to study 
the impact of turbulence and heating from both ongo- 
ing star formation, and large scale gas dynamics, on the 
Xco factor under high density conditions like the ones 
present in starbursts, mergers, and the central regions 
of galaxies. While gas density itself should have little 
direct impact on Xqq when averaging over kpc scales, 
the existence of a correlation like the SFL implies that 
higher molecular gas density is accompanied by stronger 
star formation activity. Star formation is associated with 
feedback from proto-stellar jets, stellar winds, and super- 
novae (SN) explosions, as well as heating (both photo- 
electric and by cosmic rays produced in SN), which rises 
both the temperature and turbulence of the gas. Fur- 
thermore, high density regions are also commonly asso- 
ciated with dynamically violent environments, in which 
for example, merger induced shocks and cloud-cloud col- 
lisions can induce higher levels of turbulence in the gas. 
Given the optically thick nature of the CO(l-O) transi- 
tion, any turbulence induced line broadening translates 
into a brighter I (CO) and a lower Xqo- Similarly, a 
higher gas temperature translates into a higher bright- 
ness temperature for the CO lin e. By construction, these 
effects are not recovered in the iFeldmann et al.l (j2012al ) 
simulations. 

Comparing the observed dependance of Xqo with gas 
surface density, or with CO surface brightness, to theo- 
retical numerical model predictions is challenging, and 
a series of caveats arise which must be kept in mind 
when interpreting this type of comparison. These caveats 
arise on one side from the limitations of the observa- 
tions, where the large physical spatial resolution typi- 
cally achieved in extragalactic studies 600 pc in NGC 
628 for the HERACLES beam-size) dilutes the actual 
surface brightness (or mass surface density) of GMCs in 
the beam, by some factor which depends on both the in- 
trinsic surface brightness distribution of the clouds, and 
the covering fraction of clouds across the area over which 
the emission is being integrated. In denser regions, over- 
lap of optically thick clouds can also affect the measured 
surface brightness. On the other hand, in numerical sim- 
ulations, while physical numerical quantities are read- 
ily accesible, the definition of a "cloud", and the exact 
method used to measure its surface brightness or density 
can be non-trivial. This can be especially problematic 
when the resolution of the simulation is larger, or of the 
same order of magnitude as the typical sizes of clouds. 

Narayanan et al. (2012) define Xco as the ratio be- 
tween the average column density, and the average 
CO intensity (luminosity divided by area) for each of 
their simulated galaxies, and provide a fitted formula 
to calculate Xqo from measurements of the metallicity 
{Z' = Z/Zq) and the luminosity-weighted CO intensity 
of GMCs {{I{CO)gmc)) in a galaxy: 

min[4, 6.75 x {I{CO)gmc)-°''^] x 10^° 
^co (4) 

Unfortunately, {I{CO)gmc) is not recoverable by CO 
observations at kpc scales, and also note that it is not the 



same quantity used in their definition of Xco- Instead, 
of the luminosity-weighted average of the CO intensity of 
individual clouds, we measure the average CO intensity 
over large areas of the galactic disk of NGC 628, which is 
equivalent to the average intensity going in the definition 
of Xco in the model. 

Because of averaging over large areas, with a lower 
than unity cloud filling factor, our measured I{CO) val- 
ues should be significantly lower than the CO intensities 
one would measure for individual GMCs. In order to 
use the Narayanan et al. (2012) model to make mean- 
ingful predictions from our data, we adopt a "clumping 
factor" to scale our observed average CO intensities, in 
an attempt to recover the luminosity weighted average 
intensity of individual GMC in the observed regi o ns. A 
similar approach is adopted by IKrumholz et al.l ()2009f ) 
to transform observed average II2 surface densities into 
intrinsic GMC surface densities for input to their SFL 
model. 

The choice of the clumping factor is somewhat arbi- 
trary, but it only affects the absolute value of the pre- 
dicted Xco factors, and not any relative trends observed 
across the galaxy. We use a factor of 30 to scale up 
the observed I(CO) val ues for input into Equation 10 in 
iNaravanan et al.l ()2012D . The justification for adopting 
this value is as follows. For a MW Xco factor (which is 
consistent with the average factor we measure in NGC 
628), the average CO intensity of NGC 628 inside a 7 kpc 
radius implies a mean value of Sh2 which 9 times lower 
than the typical surface density of T,h2 — 100 Mqpc"^ 
measured f or GMCs in the Milky Way and the Lo- 
cal Group (iRoman-Duval et al.l [2OIOI : iRosolowskv et al.l 
I200I IBolatto et al.ll2008ir On the other hand, compari- 
son of Equations 6 and 10 in INaravanan et all (|2012l ) im- 
plies that in their model i;//2 oc {I[CO)GMcf -^'^- There- 
fore, a clumping factor of 30 ~ gi/o.64 .^^ould imply that 
the luminosity weighted average GMC surface density in 
the inner 7 kpc of NGC 628 is similar to the typical MW 
and Local Group value. 

The predicted Xco values for different radial bins 
in NGC 628, computed using the measured metallicity 
and observed CO intensity (average between the HERA- 
CLES, BIMA, and CARMA measurements) corrected by 
this clumping fact or, are shown as blue stars in Figure [5j 
Predictions for the IFeldmann etal . (2012a) model would 
fall right on top of one of the two cyan solid lines in Fig- 
ure [5l as in their parametrization Xco is independent of 
Sg^a-nd I{CO). As stated in §5, the INaravanan et al.l 
|2012i ) model is in b etter agreement with our data than 
the model of IFeldma nn et al. (2012a). 

Figure [6] presents a comparison between the average 
Xco values for each radial bin in NGC 628 (red stars) 
and the values predicted by the INaravanan et al.l ()2012r ) 
model (blue stars). Both are presented as a function of 
I{CO) scaled by a dumpiness factor of 30, which, as 
stated above, should be roughly equivalent to the input 
{I{CO)gmc) for the model. By inspecting Figures [5] 
and [6] we can see that the model reproduces very well 
the observed trends in Xco with both metallicity and 
CO surface brightness for the inner three radial bins in 
NGC 628. The slopes of both trends are consistent, and 
the absolute values agree to Icr. 

A caveat to the above approach is that adopting a 
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Fig. 6. — ^CO a-s a function of CO surface brightness for 
the four radial bins studied in NGC 628. CO intensities have 
been scaled by a dumpiness parameter of 30. Red stars show 
our measurements and blue stars show the predictions from the 
[Narayanan et aTl 120121) model using the metallicity and CO inten- 
sity of each bin. The solid blue line shows the predictions from 
the model at constant metallicity (the characteristic metallicity at 
0.4R 25). The horiz ontal dotted line shows the canonical MW value 
from lPineda et al.l (i201fl|). 



single dumpiness parameter is equivalent to assuming 
a constant cloud filling factor over the observed area. 
A decrease in the cloud filling factor towards the outer 
disk is in fact likely, and could indeed explain both the 
flattening see n in Xgo for the N = 1 case when using 
the stacked S chruba et al.l ()2011[ ) data in the outer disk 
of NGC 628, and the discrepancy ob served between the 
data and the .Narayanan et al.l (|2012D model for the outer 
radian bin. This latter discrepancy is reduced if one as- 
sumes N = 1.5 instead of TV = 1 for the SFL. 

At this point, it is important to remind the reader 
about the many scaling factors which affect the nor- 
malization of the Xqo values. These include the zero- 
points in the flux calibration of the HERACLES CO and 
VENGA Ha maps, and the adopted Ha SFR calibration 
(we have scaled the other datasets to this reference), the 
adopted offsets used to match the metallicity scale of 
the different datasets and the models, and the adopted 
dumpiness parameter which is constrained only by de- 
manding a typical GMC surface density of — 100 
Mqpc"^. Given the inherent systematic uncertainties 
associated with all these scaling factors, it is difficult to 
draw strong conclusions from the relative offset seen be- 
tween the observed and predicted Xco values. None of 
these scaling factors affect the relative trends seen across 
different radial bins in both the data and the model. 

The main conclusion we can draw from the above 
comparison is the follo wing. Our observations are in 
good agreement with the lNaravanan et al.l (|2012[ ) model. 
If Xco follows the metallicity and CO brightness de- 
pendance proposed in this model, that is Xqo oc 
{IiCO)GMc)~°-^'^ X then, the observed metal- 

licity and CO surface brightness gradients imply that 
both quantities contribute to the formation of an Xco 
gradient across the disk of NGC 628. As stated in §5, the 
measured metallicity gradient implies Z' cx lO^*' ^'*''^, and 
a linear flt to the CO surface brightness of the inner three 
radial bins {R < 5 kpc) implies {I{CO)gmc) oc IQ-^-^'^^. 
The modd implies that Xqo « lo(0-32xo.07+o.65xo.04)i?^ 
or a gradient of Alog(Arco) = 0.05 dex kpc^^, in ex- 



cellent agreement with the measured gradient of 
Alog(Xco) = 0.06 ± 0.02 dex kpc'^ From the above 
calculation it can be seen that both the metallicity and 
the CO surface brightness (or equivalently the molecu- 
lar gas surface density) contribu t e to th e observed Xqo 
gradient. The iNaravanan et al.l (|2012[ ) prescription, in 
which A^co is inversely proportional to the metallic- 
ity, and inversely proportional to the square root of the 
molecular gas surface density, is consistent with our data. 

Our results indicate that across the disk of NGC 628, 
both the metallicity and the molecular gas surface den- 
sity (through its incidence in the temperature and tur- 
bulence of the gas) are important factors setting the CO 
to H2 conversion factor. Future analysis of a larger sub- 
set of the VENGA sample will allows us to know if this 
result is general enough to be applicable to most massive 
spiral galaxies in the local universe. 

6.2. Aco o,''^d the UV Radiation Field 

In this section we use the ionization parameter, as mea- 
sured from HII region emission line ratios in the VENGA 
data, as a proxy for the intensity of the local interstellar 
UV radiation field Suv across the disk of NGC 628. The 
goal is to study the potential effect that the UV radiation 
field has at setting the value of Aco- 

The strength of the local UV radiation field has been 
proposed as an important factor setting the CO to H2 
conversion factor. Analyzing a sample of indi vidual 
molec ular complexes in the LMC and the SMC, l lsraell 
()1997[) finds a linear correlation between the Aco and the 
ratio of the far-IR surface brightness to HI column den- 
sity {aFiR/N{HI)) which the author considers a proxy 
for the local UV radiation field per hydrogen atom (re- 
processed after absorption and reemission by dust). 

The effect of the UV radiation field on Xco is 
also discussed by iFeldmann et al.l (l2012al ) in the con- 
text of their numerical simulations and their model of 
CO emission. The authors propose that for a sin- 
gle molecular cloud, the CO abundance decreases in 
the presence of a stronger UV radiation field Ujjv 
(which they parametrize in units of the local in- 
tersteUar UV radiation field, so Uuv — S °/Smw, 

with Smw — 10^ photons cm^^ s^^ eV — 1). This de- 
pendance is stronger at low Ay (i.e. low N(H2), or low 
Z, or both), but an order of magnitude change in Uuv 
can still change the CO abundance by factors of a few 
at Ay = 1 — 3 mag. These extinction levels correspond 
to the outer envelopes CO emitting regions in GMCs, 
where CO(l-O) is starting to become optically thick. As 
the UV radiation field goes up, the CO abundance goes 
down, and Aco increases for a single cloud (see their 
Figure 2). 

When investi g ating this effect on galactic scales, 
IFeldmann et al.l (|2012aD finds that the strong depen- 
dance of Aco with the UV radiation field present in their 
model for single GMCs at fixed N{H2), is completely 
suppressed in their simulated galaxies. This is due to 
the fact that in the simulations, physical parameters 
are measured over ensembles of molecular clouds with 
a distribution in their properties, and the abundance of 
clouds with low Ay, and therefore high Aco, decreases 
strongly fo r higher values o f the UV radiation field. In 
the .Fcldmann et al.l (|2012af) simulations, increasing Uuv 
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Fig. 7. — Ionization parameter as a function of radius in NGC 
628. Black dots correspond to spaxels dominated by HII region 
emission in the VENGA data-cube. Red circles show the median 
values in 0.5 kpc wide radial bins, and red and green error bars show 
the standard deviation in each bin and the systematic uncertainty 
in the calibration respectively. 



suppresses the low density tail of the N{H2) distribution 
function, dissociating diffuse molecular clouds with low 
densities (< 10^^ cni~^) which have large Xqo values. 
Therefore, when averaging over kpc scales, the radiation 
field has little impact in the average Xqq value, which is 
dominated by the surviving higher density regions with 
N{H2) ~ 10^2 cm-2. 

A competing effect, which might be important at set- 
ting the value of Xqo, is photoelectric heating work- 
ing on large polycyclic aromatic hydrocarbon (PAH) 
molecules and small dust grains in the outer envelopes of 
molecular clouds (|Tielensll20Q5l: iGlover et aLllMot) . The 
photoelectric heating rate is directly proportion al to the 
inten sity of the interstellar UV radiation field (jTielensI 
I2005D . therefore, the gas temperature, and the CO tran- 
sition brightness temperature can be increased in the 
presence of a stronger UV radiation fiel d. By construc- 
tion, this effect cannot be recovered bv iFeldmann et all 
(p012ail as the authors impose a constant molecular gas 
temperature of lOK in their simulations. This effect 
might be present in the iNaravanan et al.l (120121 ) simu- 
lations, which follow the main heating and cooling pro- 
cesses in the gas, and are able to trace changes in the gas 
temperature, though the authors do not attempt to sep- 
arate it from dynamical effects on the brightness of the 
CO line. As mentioned above, part of the dependance 
of Xco with molecular gas surface density might indeed 
be associated with this mechanism, as higher gas density 
translates into higher star formation activity and hence, 
a stronger interstellar UV radiation field. 

In Figure [7] we present the radial distribution of the 
ionization parameter 



Sho 



(5) 



where ^'^-o is the ionizing photon flux per unit area, 
and n is the number density of hydrogen atoms. We 
have estimated q from the [OIII]A5007/[OII]A3727 line 
ratio, as measured from the VENGA data-cube of NGC 
628, following the it e rative procedures described in 
iKobulnickv fc KewlevI ()2004D . The measurements are 
limited to spaxels dominated by HII region emission 



(black dots) as in the case of the metallicity measure- 
ments described in §5.1. Red circles in Figure [7] show 
median values in at low densities (< 10^^ cm"^) 0.5 kpc 
wide radial bins, and red and green error bars show the 
standard deviation in each bin and the systematic uncer- 
tainty in the calibration respectively. 

In NGC 628, the ionization parameter falls towards 
larger radii. A linear fit implies a gradient of Alog(q) — 
0.061 ± 0.003 cm s~^ kpc~^. A drop in the ionizing pho- 
to n fiux per hydrogen atom is consistent with the results 
of lAniano et al.l (|2012[ ) who find a decreasing level of dust 
heating with radius due to a decrease in starlight inten- 
sity by modeling the Spitzer+ Herschel spatially resolved 
dust SED of NGC 628. Considering the Xco radial pro- 
file measured in §4, we find no evidence of an increase 
in Xco in regions where the UV radiation field is en- 
hanced, but we rather observe the opposite trend. The 
increase in Xc o for higher Un v expected on clouds scales 
in the model of IFeldmann et a l. (2012a ), and observed on 
molecular complexes scales bv ilsraell (jl997D . is not seen 
in kpc scales across the disk of NGC 628. This could be 
cause in part by the effects of averaging over l a rge en- 
sembles of clouds, discussed in Feldmann et al.l ()2012aD 
and described above, but this could only suppress this 
dependance. 

The fact that we observe the opposite behavior, is con- 
sistent with the photoelectric heating effect described 
above. It is difficult to decouple the impact that the 
UV radiation field might have on the gas temperature 
and hence, on the value of XcOi from surface density 
dependance discussed §6.1. All we can conclude, is that 
the observed radial distribution in the ionization param- 
eter, implies that thermal effects might be important at 
setting Xco ■ Further study will be necessary to properly 
model these thermal effects, and decouple them from dy- 
namical effects which also impact the radiative transfer 
of CO photons. 

7. CONCLUSIONS 

By inverting the SFL we obtain an independent esti- 
mate of the H2 surface density across the disk of NGC 
628 from the measured SFR surface density. Comparison 
to the observed CO intensity yields a measurement of the 
^^00(1-0) to H2 conversion factor (Xco)- By studying 
the radial profile of Xco j and its relation to other quan- 
tities like the metallicity, CO surface brightness, and the 
ionization parameter across the disk of the galaxy we 
reach the following conclusions: 

• The Xco factor increases as a function of radius 
across the disk of NGC 628. A linear fit to the data 
in the inner 7 kpc of the galaxy implies a gradient of 
Alog(Xco) = 0.06 ± 0.02 dex kpc'^ under the as- 
sumption of a = 1 SFL slope. The radial profile 
is in agreement with measurements using an inde- 
pendent technique based on dust emission model- 
ing and simultaneous fitting of the dust-to-gas ratio 
and the Xco conversion (Karin Sandstrom private 
communication) . 

• The observed Xco radial profile is independent of 
the adopted SFR tracer (Ha vs FUV-|-24/Ltm), and 
the CO emission map used (single-dish vs. inter- 
ferometer, CO(l-O) vs C0(2-l)). 
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• The observed Xqo radial profile is robust against 
changes in the adopted slope for the SFL (iV = 1 
vs 1.5). While assuming a steeper slope slightly 
steepens the observed profile, the effect is small 
compared to the changes seen in Xqq as a func- 
tion of radius in NGC 628. The observed radial 
profile is robust, and it is not strongly affected by 
potential changes in the shape of the assumed SFL 
across the galaxy. 

• The observed metallicity gradient in NGC 628 im- 
plies that regions of lower metallicity show larger 
Xco values. This is in agreement with integrated 
measurements across samples of galaxies in the 
local universe, and the predictions of theoretical 
models of CO emission in molecular gas. 

• Regions of lower CO surface brightness show higher 
Xco values. This is in agreement with theoretical 
models which predict an enhanced escape proba- 
bility of CO photons in higher density regions, due 
to an enhanced star formation activity which trans- 
lates in a broadening of the CO line due to elevated 
gas temperatures and turbulence. 

• Informe d by the theoretical model of 
iNarav anan et al.l (|2012[ ) we conclude that both 
the dependances with metallicity and H2 surface 
density contribute in roughly similar amounts to 
the formation of the observed gradient in Xco- 

• The ionization parameter shows a linear decreasing 
gradient as a function of radius, so regions where 
the intensity of the UV radiation field is higher 
show lower values of Xco- This suggest that pho- 
toelectric heating might have an impact at setting 
the brightness temperature of the CO line, and 
therefore the value of the CO to H2 conversion fac- 
tor. 

• When comparing with theoretical models, our ob- 
se rvations agree ve r y we ll with the predictions 
of INaravanan et aD ()2012). While our d a ta also 
agrees with the model of iFeldmann et al.l ()2012a[ l 
in the slope of the Xco-metallicity relation, and 
the lack of dependance of Xco with the local ion- 
izing field on kpc scales, this model predicts Xco 
values which are typically a factor of ~ 4 higher 
than the observed ones. 

• Given the observed Xco radial profile in NGC 628, 
we conclude that using a single MW Xco factor to 
estimate the total H2 mass would imply an overes- 
timation of 20% in this quantity. While this is a 
relatively small effect, our results imply that much 
larger systematic deviations can occur when using 
a canonical value in specific regions within galaxies, 
or when the local conditions of the ISM in the re- 
gions of interest differ significantly from the canon- 
ical values. 

A series of assumptions have been made in order to 
reach these conclusions, and it is important to keep these 



in mind while interpreting our results. The uncertainty 
in the flux calibration and sensitivity to emission on dif- 
ferent spatial scales in the three CO datasets used, the 
use of a constant non radially dependent C0(2-l) to 
CO(l-O) ratio to scale the HERACLES data, the un- 
certainty in extinction correction (in the case of Ha) and 
the SFR calibrations used for Ha and FUV plus 24^m, 
and the assumption that a single molecular SFL holds 
across the disk of the galaxy, all add up to make the sys- 
tematic uncertainty in the absolute value of Xco values 
large. An important assumption made is that all these 
uncertain scaling factors do not have a strong radial de- 
pendance, therefore allowing for the relative radial trend 
in Xco to be measured. This assumption is supported 
by the fact that the six combinations of datasets used 
in this study, as well as independent measurements from 
dust SED modelling (Karin Sandstrom private commu- 
nication), show similar trend of an increasing Xco with 
radius. 

In the future, we expect to extend this study to a 
larger sample of galaxies in the VENGA survey, in 
order to confirm if the observed trends are a common 
feature among massive spirals in the local universe. In 
the near future, ALMA will allow us to study molecular 
complexes in nearby galaxies like NGC 628, in the same 
level of detail that we can currently achieve in the MW. 
In the context of this study, the commissioning of the 
MUSE IFU in the Very Large Telescope (VLT), wih 
provide an instrument with excellent spatial resolution 
over a large field-of-view, which could be the perfect 
complement for ALMA observations of nearby galaxies. 
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